Radiation-reaction in classical offshell 
elect ro dynamics : 
I. The above mass- shell case 

I. Aharonovich 
L. P. Horwitz 

January 26, 2013 

Abstract 

Offshell electrodynamics based on a manifestly covariant off-shell rela- 
tivistic dynamics of Stueckelberg, Horwitz and Piron, is five-dimensional. 
In this paper, we study the problem of radiation reaction of a particle in 
motion in this framework. 

In particular, the case of above-mass-shell is studied in detail, where 
the renormalization of the Lorentz force leads to a system of non-linear 
differential equations for 3 Lorentz scalars. 

The system is then solved numerically, where it is shown that the mass- 
shell deviation scalar e either smoothly falls down to (this result provides 
a mechanism for the mass stability of the off-shell theory), or strongly di- 
verges under more extreme conditions. In both cases, no runaway motion 
is observed. Stability analysis indicates that the system seems to have 
chaotic behavior. 

It is also shown that, although a motion under which the mass-shell 
deviation e is constant but not-zero, is indeed possible, but, it is unstable, 
and eventually it either decays to or diverges. 

1 Introduction 

Classical 5D electrodynamics arises as a £7(1) gauge of the relativistic quan- 
tum mechanical Stueckelberg-Schrodinger equation (shown below in (3)), sim- 
ilar to the construction of Maxwell fields from the U(l) gauge of the classical 
Schrodinger equation [14, 19, 20, 27]. 

The fields depend on the spacetime points as well as the invariant evo- 
lution parameter r, and thus, are defined on a five-dimensional manifold. In 
previous papers, we have studied the configuration of these fields associated 
with a uniformly moving source [1] and uniformly accelerating one [2]. 

In this paper, we continue the investigation of these fields and examine the 
problem of radiation-reaction of a point-particle interacting with its own field, 
using the explicit r-retarded Green- Function first presented in [2]. 



In this relativistic framework, the mass of the particle m 2 — —p^ (or 
its gauge invariant generalization in the presence of electromagnetic fields) is a 
dynamical variable, whereas M, which enters into the Hamiltonian formulation 
of the theory (see below equation (2)), the so-called Galilean target mass[-i] 
(defined as the on-shell mass), is the non-relativistic limit of m. One then 
defines a dimensionless mass-shell deviation scalar 

m 2 - M 2 



which provides a measure of how far the particle is offshcll. 

As the fields are essentially 5D, the Huygens' principle no longer applies (cf. 
[9, 17, 18, 21]), and the r- Green- Function given in ref. [2] has support on the 
entire history of the particle, within its local 5D past lightcone. Therefore, the 
Lorentz force for the self-interacting particle with its own field, depends on the 
particle's own entire history within its past lightcone 1 . 

It is shown in section 3 that the radiation-reaction force differs substantially 
for the three possible cases of e, i.e., the above-mass-shell case e > 0, which is 
the main topic of this paper, the on-mass-shell case e = 0, and the below-mass- 
shell case e < 0. This is mainly due to the 0(4, 1) choice of symmetry for the 
homogeneous fields 2 . 

For the above-mass-shell case e > 0, it is shown that the force has a renor- 
malizable simple-pole singularity. After renormalizing this pole, the residue is 
a local vector equation with up to 4 th order terms in x^. After finding an ex- 
pression for a;' M which depends only on lower order terms ± M ,i M and we 
are able to construct a set of 6 coupled scalar equations, with three scalars and 
their derivatives as dynamical variables, where e, the mass-shell-deviation, is 
the leading scalar. 

Apart from the constant e case, the analysis to be shown below is numerical 
(see section 4), performed using the high-precision math library MPFR [8]. As 
initial condition, e begins with a non-zero positive value, as the particle is above 
its mass-shell; it either decreases smoothly towards e — > 0, or diverges to e — > oo. 
For the e — > case, the particle relaxes to an asymptotically on-shell uniform- 
motion. In the e — > oo case, the acceleration d 2 x^ / At 2 increases indefinitely. 
It is shown, however, that the 3 velocity dx l / dt converges strongly to a finite 
value, such that in any particular Lorentz frame, the particle is asymptotically 
inertial. 

Therefore, in both cases, there is no runaway motion in x(t). Furthermore, 
the state of motion for e — const ^ is shown to be unstable. 

For the noninteracting Stueckelberg Hamiltonian, there appears no a priori 
reason for the particle to have a particular mass. The above-shell analysis of the 
particle in self-interaction shows that for a wide range of conditions converges to 

1 In this paper we do not consider external fields. 

2 In [1], for example, the choice of initial mass shell for the uniformly moving source was 
shown to lead to substantially different pre-Maxwcll fields, even though their zero-mode cor- 
responds to the same Maxwell fields. 
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the on-shell value. This striking result is an important indication of the stability 
and physical applicability of the Stueckelberg dynamics. 

The remainder of this paper is organized as follows: In section 2 an overview 
is provided of the manifestly covariant relativistic dynamics of Stueckelberg, 
and the corresponding offshell-electrodynamics. In section 3 the Lorentz force 
is formulated using explicit r- Green- Functions [2], and it is shown how the 
sign of the mass-shell-deviation e leads to substantially different results. In 
this paper, the above-mass-shell case is analyzed in detail to provide equations 
of motion (the above mass-shell condition is preserved by the motion). We 
shall deal with the below-mass-shell problem in a later paper. In section 4, the 
equations of motion are integrated, and the various results are shown. In section 
5, conclusions of the analysis and simulation are given, with prospects of further 
investigation. Furthermore, the regularization and renormalization procedures 
of Gel'fand [10] are given for reference. 



2 Fundamentals 

2.1 Relativistic dynamics 

Stueckelberg postulated a manifestly covariant Hamiltonian-type evolution in 
both classical and quantum mechanics [29, 30], by introducing the Hamiltonian- 
type evolution function for a free particle 3 , 

K = j^fcp" (2) 

Phase space is 8D and consists of = (t, x) and = (E,p), developing with a 
Lorentz invariant parameter t. Particles in this framework trace out worldlines 
in spacetime, and are therefore, also denoted as events. 

The corresponding quantum equation of motion is the Stueckelberg-Schrodinger 
equation (we take h — 1) 

i^ T {x) = K^ T {x) (3) 

The framework was extended to a many-body system by Horwitz and Piron 
[12], where r was given a physical significance as a universal historical time, 
essentially that postulated by Newton. This interpretation was necessary in 
order to write equations for the many-body system. The evolution function K 
can then take a more general form: 

N 1 

K = WTfV^PnPn + V(X1,X 2 , ■■■ 1 X N ) (4) 

71=1 " 



3 Throughout the text, we shall work with standard relativistic notation where fi,u £ 
{0, 1, 2, 3}, and the Lorentz metric signature has the form (—1, 1, 1, 1), and we take units for 
which c = 1 
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V(x\,X2, ■ ••) = V(x) is a Lorentz scalar function of spacetime positions of all 
the particles at the (universal) time r, which establishes a dynamical correlation 
between them. The classical equations of motion are similar to those of the non- 
relativistic Hamilton equations 

dp ntl M n dx% ox% 

In the usual application of special-relativity (SR, cf. [25, 28]), the energy- 
momentum is constrained to a mass-shell defined as: 

p"P M =P 2 -E 2 = -M 2 (6) 

where M is a given fixed quantity, a property of the particle. In the Stueckelberg 
formulation, however, the particle's mass is a dynamical property which may 
depend on r: 

p% = -m 2 (7) 
The relation between r and proper time s is given, according to (5), by 

1 777 2 

ds 2 = - dx^ dx = -x^x dr 2 = dr 2 = - dr 2 (&) 

Thus, the proper time interval ds n and the universal time interval dr, are related 
through the ratio between the dynamical Lorentz invariant mass m n , and the 
so called Galilean target mass M n . If V(x) — > asymptotically for large r, then 
each 77i 2 becomes a constant, since dp£/ dr — > 0, we find 



K = - vtVn a = — ^ — — = const 



Since this asymptotic value is what we expect to measure in an experiment, we 
may take M n to be this asymptotic value. Is has been an open question in the 
framework of this theory as to why a given particle appears (even approximately) 
to have the same value of m in many circumstances, e.g., after many scatterings. 
Possible answers lie in the properties of a particle in interaction with other 
particles 4 , or a relaxation to a minimal free-energy. Some insight has been 
gained in the framework of relativistic mechanics [13] in the discovery of a high- 
temperature phase transition for an offshell particle to become confined to a 
particular mass shell [5]. In this work, a new, renormalized mass M reg emerges, 
related to M by some very large scale, and a new, very natural and intrinsic 
mechanism is found, in which the self-interaction of a relativistic charged particle 
can cause a relaxation of the particle mass to its renormalized mass. A possible 
relation between these results will be studied elsewhere. 

4 T. Jordan, personal communication. 
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2.2 Off-Shell Electrodynamics 

What we shall call "pre-Maxwell" off-shell electrodynamics is constructed in a 
similar fashion to the formal construction of standard Maxwell electrodynamics 
from the Schrodinger equation [27] (see also [32, 33]). 
Under the local gauge transformation 

%(x) =e- teoX< - x ^^ T (x) (9) 

5 compensation fields a a (x,r) (a <E {0,1,2,3,5}) are implied by (2) and (3), 
such that with the transformation 

a ' a ( J) - T ) = «« {x, t) - d a x(x, t) 

in the following modified Stueckelberg-Schrodinger equation remains form in- 
variant 



i— + e a 5 (x,T) 

OT 



*t(x) = ^ [(p" - eoafXpv - eo^)] * r (x) (10) 



under the transformation (9). 

The resulting K evolution function is then of the same form as for the usual 
U(l) gauge compensation argument for the non-relativistic Schrodinger equation 
(where the fourth field a arises from the requirement for the gauge invariance 
for the time dependent Schrodinger-equation) , here given by: 

1 2 
K = \P ~ e oa(x, t)] -e a J (x,r) (11) 

(where we have used the shorthand notation of M ), and the correspond- 

ing Hamilton equations are 

*.(r) = g = ^-eo^] (12) 

r{r) =-^ = Jj(P- e oa(x, t)) v &>a v (x, r) + eo&>a*(x, r) (13) 

Here, eo is proportional to the Maxwell charge e through a dimensional con- 
stant, which is discussed below. Second order equations of motion for ^(t), a 
generalization of the usual Lorentz force, follow from the Hamilton equations 
(12) and (13) 

MF = eox v f" v + e r 5 (14) 

where for a, /3 = 0, 1, 2, 3, 5 the antisymmetric tensor 

f af} = d a a P - d p a a (15) 

is the (gauge invariant) 5D field tensor. Moreover, second order wave equation 
for the fields f a & can be derived from a Lagrangian density 

Se = -\f a pf aP -e a a j a (16) 
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which produces, by the usual variational methods, the wave equation 

Xd a f a = e Q f (17) 

A is a dimensional constant, which will be shown below to have dimensions of 
length. The sources j^(x, r) depend both on spacetime and on r, and obey the 
continuity equation 

d a j a = drf" + d rP = (18) 

where j 5 = p is a Lorentz invariant spacetime density of events. This equation 
follows from (10) for 

p T [x) = %(x)^ T {x) 

=~ [KM (id" ~ ^(x, r)) * T (z) + c.c] 
as we discuss below, and also from the classical argument given below. 
2.2.1 Currents of point events 

Jackson [16] showed that a conserved current for a moving point charge can be 
derived in a covariant way by defining the current as 

J fl (x) = e dsz^(s)5 4 [x~ z(s)} (19) 



In this case, s is the proper time, and z^(s) the world-line of the point charge 
(for free motion, s may coincide with t), and i M (s) = — z M (s). Then, 



r+oo j 

da J M = -e / ds — 5 4 \x - z(s)} = -e lim 5 4 [x - z(s) 



(20) 



which vanishes if z^(s) (or, for example, just the time component z°(s)) becomes 
infinite for s — > ±oo, and the observation point x^ is restricted to a bounded 
region of spacetime, e.g., the laboratory. We therefore, with Jackson (see also 
Stueckelberg [29]), identify J' 1 as the Maxwell current. We see that this current 
is a functional on the world line, and the usual notion of a "particle" associated 
with a conserved 4-current (and therefore a conserved charge corresponding to 
the space integral of the fourth component), corresponds to this functional on 
the world line. 

If we identify S 4 [x — z(s)] with a density p(x, s) and the local (in r) current 
z fi (s)5 4 [x — z(s)] with a local current j^(x, s) 

p(x,s) = 5 4 [x - z{s)} f(x,s) = z^(s)S 4 [x-z(s)} (21) 
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then the relation 



^-J 4 [x - z(s)] = -z^{s)d^S 4 [x - z(s)] 

used in the above demonstration in fact corresponds to the conservation law 
(reverting to the more general parameter r in place of the proper time s) (18) 

d li f{x,T) + d T p{x,T) = {) (22) 

What we call the pre-Maxwell current of a point event is then defined as 

j a (x,T) =z a (r)5 i [x-z(T)\ (23) 

where j 5 (x,r) = p{x,r) and i 5 (r) = 1 (since z 5 (t) = t). Note that j a does 
not have 0(4, 1) tensor properties, since the factor i M in cannot be obtained 
from j 5 by a linear transformation (this also occurs in the 4D theory where 
the current is constructed from a 5 3 (x) for the J° component and the three- 
vector components J 1 with velocity times this S 3 (x) 5 ), and therefore, the fields 
associated with this current are not tensors under 0(4,1); the components 
are, however, tensors under 0(3, 1), the Lorentz group. 

Integrating (17) over r, we recover the standard Maxwell equations for 
Maxwell fields defined by 

A fl (x) = J a»(x,T)&T (24) 

We therefore call the fields a M (x,T) pre-Maxwell fields. We see this by noting 
(16) 

/-foo Z'+OO 
[d^(x,r) +d 5 F 5 (x,T)] dr = A / d a r a (x,T)dr 
-co J — oo 

/+oo 
jp(x,r) dr 
-OO 

= e a h»(x) (25) 

where the r integral of d^f^ is assumed to vanish. This can be accounted 
for by considering the field as a wave-packet that propagates to infinity in 
space time for r — > oo, or, in terms of its Fourier transform to its offshcll mass 
distribution, 

/+oo 
e^a^r) dr, (26) 
-co 



5 The covariant Maxwell theory assumes J^(x) to be a given 4- vector (not based on the 
Jackson type construction [16]) and therefore, the theory is of course, Lorentz covariant. We 
have no a priori reason however, to assume 0(4, 1) covariance of the theory. 
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as an application of the Riemann-Lebesgue lemma. Eq. (24) then implies 
that the Maxwell potentials and fields correspond to the zero-mode of the pre- 
Maxwell offshcll fields, with respect to r, i.e., 

A"(x)=a»(x,K)\ K = (27) 
The relation (24) then follows. Moreover, since 

d v F>^ = ^J^(x), (28) 
Ae 

the constants A and eo are related by 

A = 7 e= X (29) 

where e is the standard Maxwell charge. 

For the quantum theory, a real positive definite density function p T {x) can 
be derived from the Stueckelberg-Schrodinger equation (3) 

P(x,t) = \V t \ 2 = %(x)* T (x) (30) 

which goes over to p{x,r) = J 4 [x — z(t)\ in the classical (non-quantum) limit. 
It follows from the Stueckelberg equation (10) that the continuity equation (22) 
is then satisfied for the gauge invariant current 

f{x, r) = ~ [KWW ~ e o^(x, t))* t (x) + c.c] (31) 
Using (16), applicable in the quantum theory as well [20, 27], and (24), we find 

/■+0O 

J"(x) =e / j"(x,T) dr (32) 



2.2.2 The wave equation 

From equations (17) and (15) one can derive the wave equation for the potentials 

a a (x, t): 

\dfsd p a a - A(9 Q (a /3 a^) = e j a (33) 

Under the generalized Lorentz gauge dpa? — 0, the wave equation takes the 
simpler form 



\dpd p a a = A 



B 2 n a 



e Q j a (x,r) (34) 



where the 5 th diagonal metric component can take either sign cr 5 = ±1, corre- 
sponding to 0(4, 1) and 0(3, 2) symmetries of the homogeneous field equations, 
respectively. 
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Integrating (34) with respect to r, and assuming, as above, that lim a a (x, r) 
we obtain 



+ OG 



dr 



2 a a + cr 5 



d 2 a a 
Ih 2 



eo 



r{ x ) 



With (24), we obtain 



e 



(35) 



i.e. 



a 2 A^(x) = J»{x) (36) 

(for (jt = 0,1,2,3) 

Therefore, the Maxwell electrodynamics is properly contained in the 5D 
electromagnetism. 



2.2.3 A note about units 

In natural units (H = c = 1), the Maxwell potentials have units of 1/L. 
Therefore, the pre-Maxwell offshell potentials a a have units of 1/L 2 , and in 
order to maintain the action integral 

/+oo 
JS? dr d 4 x (37) 
-oo 

dimensionless, the coefficient A in (16) must have units of L, forcing eo to have 
units of L as well (e is dimensionless). 



2.3 Solutions of the wave equation 

In our previous works [1, 2], we have discussed the Green- functions associ- 
ated with the 5D wave-equation (34). In particular, using the ultrahyperbolic 
Riemann-Liouville integro-differential operator of Nozaki [22], the r-retarded 
form was shown in to be (for the 0(4, 1) case, which we shall concentrate on 
here, a 5 = +1): 

, , 6M 9(-x 2 -t 2 ) 
^ T) = -4 ^-V 2 (38) 



where: 



Ju I Ju ft Ju I Ju (y Ju 



Thus, the fields of a point source can be directly integrated to yield: 

a a (x,T) = e J dV dT'g(x-x',T~T')j a {x',T') (39) 
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and applying it to a point particle represented by (21). The potentials a a (x, r) 
are then given by: 

/+oo 
At' g(x - z{t'),t ~ t') z a {r') 
-OO 

= -^nf o[-(*-^)Y-(t-t>Y] , dr , m 

4tt 2 ./_«, [-(x - z(t')Y - (r - T') 2 ] 3 / 2 1 ' [W) 

where the 1Z(. . .) indicates that the integration (40) is regularized at the points 
where the denominator vanishes, which, by the 9{. . .) function in the numera- 
tor, are also the bounds of the integration. We shall explain this regularization 
procedure in detail in the following. 

The r-retarded Green-Function (38) leads to r-retarded solutions of the 
wave-equation. However, their zero-mode, identified with the Maxwell field 
(27), is not t-retarded, since the r-integration of the fields (40) are, essentially, 
symmetric under time i-reversal, and results in the average of i-retarded and 
i-advanced potentials. This average does not contribute to the Maxwell radia- 
tion field [7], and is singular in the zero-mode. It is, however, propagating on 
the pre-Maxwell level and well defined by (40). 

Equation (40) above is the foundation on which the Lorentz force of back- 
reacting particle is evaluated in the next section. 

3 Lorentz force 

The potentials a a (x,r) (40) can be used to find the corresponding fields f a p, 
and with them, the Lorentz force (14). 

We then find the radiation-reaction equation to be 

Mx!>{t) = (-l) 2 ^li a (r)n £j(R(ry))^^dr' (41) 

where 

R(t, t') = -(x a (r) - x a (T'){x a {T) - x a (r r ) = 

= -(x,(r) x„(t'))(x»(t) - x^t')) - (r - r') 2 (42) 

The integration is bounded by regions where i?(r, r') is positive within the 
range — oo < r' < r, i.e., within the past 5D light cone of the particle at x a (r). 
The numerator function h a ^(r, t') is given by: 

h af3 (T,r') = x a (T>)-\-R(T,T>) - xV{t')-\-R{t,t>) (43) 
ax/3(T) ox a (T) 

Defining h = r — r' , we expand R(t, t') — R(t, t — h) around h = 0: 

h 4 h 4 

R(t, t') = -h 2 x a x a + h 3 x a x a - —x a x a - 2—x a x a + 0(h 5 ) (44) 
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Similarly, expanding h a ^(r, t — h) around h = 0: 

^(r, r') = ft 2 (± a a^ - ±<V) + 0{h 3 ) (45) 

Let us define the leading term in R(t,t') as e: 

e(r) = -x a x a = -x^ - 1 (46) 

where we have used the definition of x 5 (t) = r and x 5 (t) = 1. To see how this 
coincides with the definition given in the introduction (1), let us write: 



M 2 M 2 



and thus, clearly: 



m 2 m 2 — M 2 

e = — x^x -! = — -! = 



M 2 M 2 



The function J2(r, t') is the generalized Lorentzian distance function (interval 
on the 4 + 1 space) between 2 points in the particle's history, namely, its present 
location at a; a (r) and its previous location at x a (r') at the time r'. The distance 
converges to R(t, t') — » as t' — > t, or as /i = r — r' — > + . However, the sign of 
e at t determines whether the immediate past is inside the domain of integration: 

• We call the case of e(r) > the above mass-shell case. In this case, the 
integration runs up to the point r = r' itself, where there is a strong 
singularity. From the leading terms in R(t, t — h) and h a @ (t, t — h), we 
see that the singularity is of the order of 2-5/2 — 2 = 3. 

• e(r) < is the below mass-shell case, where the integration is bounded 
by some T\ < r, where R(t,ti) = and R(t,t — h) < for all < h < 
r — T\. I.e., in such a case, the particle is 5D spacelike with respect to its 
immediate past, up to t\, and therefore, that part of the history does not 
influence, or back react, with the particle at its present position. 

The interaction is only with parts of the particle's history which are inside 
the 5D past cone with respect to its present 5D position x a (r). 

In the neighborhood of these crossing points, denoted by n, the expansion 
R(t, rj) takes the form: 

R( T , n-h) = hR(n) + 0(h 2 ) (47) 

I.e., for the general case, R{t{) ^ 0, i.e., where the motion in and out 
of the past light-cone is well defined, the singularity due to the vanishing 
denominator is of order 5/2, namely, non-integral order, which enables 
regularization to remove it. 
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• e(t) = 0, is the on mass-shell case, m 2 — M 2 . In this case, higher order 
terms determine whether the orbit in the immediate past t' — > r is inside 
the past 5D cone or not. Locally, at t' — r, the particle's orbit is tangential 
to the 5D cone. If this is persistent, i.e., if the on-mass-shell motion existed 
over some measurable interval in the particle's history, and the immediate 
past is inside the lightcone, then the singularity of the denominator is even 
stronger, at least of order 4-5/2—2 = 8. 

The various possible trajectories are shown in figure 1, where, in this case, the 
on-shell motion is actually outside the lightcone, then becoming tangential at 
t' = r, or h = 0, at the vertex of the cone at x^t). 

Henceforth, we shall study the above mass-shell case, and thus, assume 
e(r) > for all the range of relevant t. The other cases will be studied in 
a succeeding publication. 

Let us begin by breaking up the integration (41) as follows: 



/ 



8tt 2 



r~ h a/1 (T r'l ( T h afJl (r i 

(48) 



Let us focus on the second term, the integration range r' G (r — A,r). A 
is chosen such that R{t,t') ^ for all r — A < t' < r. For most practical 
purposes, we shall take A = 1. In any case, its non-zero value can be absorbed 
in the renormalized mass, to be defined below. 

In the range (r — A,r), we can therefore expand both numerator and de- 
nominator. Since e ^ 0, and both h a ^{r, t') and R(t,t') are 0(h 2 ), we can 
define the following functions: 

h a »{T,T') EE h 2 X k a "(T,T') (49) 

R(T,r') = h 2 xT(r,r') (50) 
We then find, for the second term in (48): 



_ 3e 2 ± a (r) r h^(T,r') ,_ 
_ 3e§i a (r)^ f A h 2 xk^( T ,T>) dh 



87T 2 J (/,2)5/2 T 5/2( 7 

_ 3egi; a (T) ^ f A fc^(r,r') dh _ 3e 2 i a (r) „ f A ^{r,h) 
-^^ n ] TV%7)h? = ^^ Jo ~h^ dh (H) 

where we have changed integration variable to r' — > h, and defined the function 

^(T, h) EE k a ^ l (T,T-h)/T 5 / 2 (T,T-h). 

If we consider <fi afl (t, h) as a test function, and 1 /h 3 , as a generalized function, 
we can apply Gel'fand regularization [10]. Gel'fand showed that the generalized 
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Figure 1: The past light cone, where the vertex in the origin is the present point 
x a (r) (only t,r and r axes are plotted). The 5D light cone interior is given by 
(£(t)— t(r')) 2 > (T~T') 2 + (r(T') — r(T)) 2 , where the axis for h = t — t' is drawn. 
Retardation is when h > 0, bounded by the hatched surface at ft. = 0. 3 possible 
trajectories are shown: inside the past lightcone, plotted in green, corresponding 
to the e(r) > case, where e(r) is essentially the slope i(r) 2 — r 2 — 1. The case 
of e(t) < is shown in red and e{t) — in blue. 
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function h+ 

h\ = { kX h> ° (52) 

|0 h<0 

has simple poles for all negative integers A = — 1, — 2, . . .. The residue at these 
poles is given by 

A *_„// = (53) 
(n 1). 

Returning to the full equation of motion 

3egi Q (r) r r~ A ^(r,r') , , , f A k^(r,r') dh 



IF 



8tt 2 



71 / ' ' dr' 

-co 4 /2 (r,r') Jo T^(t,t') 



(54) 



If the first integral (54) is regularizable, then it is finite. The second integral, 
though non-regularizable, is renormalizable. The renormalization procedure 
works by dividing both sides by Y{— 3 + e), and taking the limit of e — » 0. 
However, in order to ensure that the equation of motion still maintains a mass 
term, we shall renormalize the mass as well. Let us then define 



M 

Mrcg = }™ fP3Te) (55) 



If M is finite, then M rog — !• 0. We therefore assume that m — > oo such that 
M xes =/= 0, as is often the case in renormalization in quantum field theory, i.e., 
in the abovc-mass-shcll case, it is the renormalized mass that enters into the 
dynamics. 

Recall that e depends on the ratio of m/M; since e is finite, and determined 
by the dynamics, it follows that m is also renormalized. This implies further that 
the quantity p^ — eoa^ is essentially renormalized in its scale, but i M remains 
covariantly and finitely determined by the dynamics. In the stable situations 
that we describe below, when e — > 0, m reg — > M rcg , and it is this value of the 
mass parameter that may be identified with the Galilean target mass. 

On the other hand, the first integral (54), presumably regularizable, vanishes 
under the renormalization, and thus, we are left with the renormalized mass and 
the residue of the second integral, given by: 



„ _ 3e 2 x a (r) d 2 ^(r,fr) 



M r 



(56) 

h=0 



We therefore, need to expand the components of <p a ^ up to h 2 . Expansion 
reveals: 

fc^(r, r-h) = bf + hbf + ^h 2 bf + 0(h 3 ) (57) 
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where: 





= x a x p - x p x a 


bf 


4 

~ 3 


x a x - x^x 01 


b ap 


CO 1 CM 





&*3p - afx" 



(58) 



Similarly, for T(r, r — h) we find: 



T(r, T-h)=r Q + h ri + -h 2 r 2 + 0(h 3 ) 



(59) 



and its expansion terms: 



ro = — x a x = e 

■ -a 1 

ri = x a x =--t 

3. .... : 



'2 



_ *Xj r\ tX^ _ *Xjr\tjC 

2 2 



(60) 



Using the Axiom computer algebra system [24], the second derivative of 



4> a ' (t, h) at /i = is found to be 

d 2 Q ' 3 (T, ft) _ -lOfoo' 3 r r 2 + 356o"r? - 206? fl r n + 4bfr% 



dh 2 



4r, 



9/2 



(61) 



Defining 



M rGg ■ 8tt 2 
3e 2 



(62) 



we finally obtain the renormalized Lorentz force equation: 



Dx^ 



/ 45 35 \ 

i M ( bees H eep i 3 + 2e 2 x a x OL ) 

V 6 8 / 



/ 35 \ 

+ x" ( - 4e 2 e - 3e 2 p + y^ 2 ) - Se^P 1 + 2e 3 iF 



(63) 



where we have defined an additional scalar p given by: 



(64) 



where we note that for any derivative higher than unity, the fifth component of 
x 5 is annihilated, as x 5 (r) = r and x 5 (r) = 1 . 
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This Lorentz force equation is of 4 th -order. It is more natural to express the 
higher order part in terms of lower order. From equation (63), we find that 



i M (beee 



........ 45 . 35 

x^ I beee + —eep — — s 
6 8 



55" ( - 4e 2 e - 3e 2 p + ^ee 2 - D x 8e 9 / 2 ) - 6e 2 ex» 



n (65) 



where 

M"„ = e^ v + x"x v (66) 
One can eliminate the factor „ by contracting with x^: 



(67) 



Thus: 

x l^^t u X — X i/ X 

2 



■ -u(r ■■■ 45 • 35 -3\ 

x^x* [beee + — eep - — e ° J 

+ i M x"( - 4e 2 e - 3e 2 p + yee 2 - £> x 8e 9/2 ) - 6e 2 ei M x M 

^ ••• 45 • 35 • 

(-£ - 1) I 5eee + —eep — — e' 
V 6 8 



ie( - 4e 2 e - 3e 2 p + yee 2 - D x 8e 9 / 2 ) - 6e 2 ei M ^ 



(68) 



We then find a simpler equation of motion for x v : 



.„ ./35 . 



4£>e 7/; 



35 \ 
4ee + 3ep - —e 2 + 8De 7/2 J + Qeex^ 



(69) 



3.1 Scalar set of equations 

Equation (69) is a vector equation of 4 th order, up to 3 rd derivative of x^ . 
Therefore, the of this system is 12D, if one excludes external fields, comprising 
the dynamical variables ±",55" and x M . If external fields are also present, then 
x 1 - 1 is also a part of configuration-space, increasing the dimensionality of the 
system to 16D. 

Reduction of phase space is possible by suitable contractions of (69). 
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Ki 


2i 
~ e* 


r35 , 


- ADe 7 ' 2 


K 2 


2 


4ee + 




K 3 


= 12 £ - 

e 





For this purpose, we define yet an additional scalar, 77, 

77 - (70) 

Using the definitions of e, p, 77 and their derivatives, and using the following 
relations 

x li 'x> i = -ig - ~p x M x M = - 77 = ^77 

we eventually find the following set of equations: 

e = -3p + 2(e + + £#2 + (e + 2p)if 3 

p = 277 - eXi + 2p J fY 2 + pX 3 (71) 

77 = -(e + 2p)isTi + pX 2 + 277^3 

where we have defined the following auxiliary functions, for clarity: 



(72) 



Equations (71) indicate that the reduced phase space includes the scalars 
{e, i, e, p, p, 77}, all of them containing up to 3 rd order derivatives in x^. Numer- 
ical study of these equations is given in the following section 4. 

The Ki functions (72) can be interpreted as potential scalars, as they prac- 
tically form the potential landscape of the system. 

Before that, we can still study two simple solutions analytically, both of 
them, for constant mass-shell deviation e. 

3.2 The motions with constant mass-shell deviation 

In this section, we are seeking solutions such that e remains constant. Thus, 
this immediately leads to: 

i = £ = £ = 

In this case, only K 2 = 6p/e + 16De 3 ^ 2 remains non-zero, and the only non- 
trivial equation is the one for p: 

p' = 277 + 2p(6^ + 16L>£ 3/2 

However, from the equation for e — 0, we immediately find p — 0, and therefore, 
we are left with: 

£ = £=£ = /3 = 77 = 

r, = - p (*E + l6De^ (73) 
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In this case, all e,p and rj remain constant. 

Clearly, the simplest case would be to set p = r\ = 0, which is the case of 
uniform motion. This is understandable, as uniform motion should not cause 
radiation reaction. 

The other case would be a non-zero p, 77, which means that even for acceler- 
ated particle, a particular orbit exists for which e remains constant. 

Both cases, as seen below, are actually unstable, and the case of non-zero 
p, 7/ is unstable numerically as well, eventually leading the particle either back 
to mass-shell e — > or to infinite mass e — > 00. 

4 Numerical simulation and results 

In this section, numerical simulation of the set of equations (71) is shown. The 
scalar system comprises the 6 scalars e, i, e, p, p and 77. The scalar system offers a 
reduced phase-space description of the system, and its divergence is completely 
unrelated to a particular Lorentz system. However it lacks the coordinate repre- 
sentation offered in the original vector system given in (69). Therefore, when it 
aids clarity, this vector equation set was also used to produce numerical results. 

4.1 Numerical methods 

The system was studied both in the scalar set of equations (71) and the original 
vector one (69). The simulation was written in CH — h, employing a standard RK4 
integrator 6 taken from GSL [ : ] , and rewritten to use high-precision arithmetic 
provided by MPFR [ ] and MPFRCPP [ ] as its higher level interface. Results 
are plotted mostly with matplotlib [15], and some with Visit [6]. 

All the numerical results of this paper were obtained with free and open- 
source software. 

4.2 Results 

Let us begin with the simpler cases, where e smoothly goes to 0. In figure 2, 
3 plots of e(r) for 3 initial conditions for £0 are shown. The other variables, 
namely £o,po,po and 770 all begin at 0. Even though e tends to increase e, 
eventually, s, with all derivatives, smoothly go to zero. 

In figure 3, the entire set of scalar variables showing the entire behavior of 
the system, as it smoothly falls down towards the mass shell and asymptotically 
into unifom motion. 

We note that the asymptotic limit of e = cannot be handled in this set, as 
the point e = is not defined for this system. Therefore, in order to handle the 
transition between above-mass-shell, on-mass-shell and below-mass-shell behav- 
ior, one has to develop and solve (41) for the other cases as well. This would be 
studied in a succeeding publication. 

e Runge-Kutta method of 4 th order with 5 th order error estimatation. 
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Figure 2: An example of 3 solutions of (71), where the mass-shell deviation 
scalar e(r) starting from non-zero value and eventually falling smoothly towards 
the on-shell e — > 0, asymptotically approaching uniform motion. 
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Figure 3: A converging e — > example, where the full set of scalar variables is 
shown, with initial conditions e = 0.5, Sq = 0.1, e'o = po = Po = f]o = 0. 
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Figure 4: A converging e — > example, using the vector system of equations 
(69). Initial conditions are given in vector form for Xq,£q, and where Xq = 0. 
The entire motion is in the t — x plane, and clearly, both x and t converge 
towards uniform motion. 
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In figure 4, a similar system is solved using the vector system of equations 
(69), with initial conditions in x^,x^ and x^. It is evident that as e — > 0, both 
t(r) and x(t) approach a constant value, and as a consequence, v x = x/t — 
dx/ dt goes to a finite limit as well (less than one). 

In figure 5, a strongly divergent solution of the scalar system is shown. In 
this case, p is negative, whereas in ordinary SR dynamics, we find, for an on- 
mass-shell particle, that x^x^ = a 2 f > 0, where a\{ is the acceleration observed 
in a co-moving frame. In this particular case shown in 5, the initial po was 
already negative, and therefore, cannot be associated with a simple 3-vector 
acceleration in a co-moving frame. It is possible, however, to arrive at a 
diverging case, even when initially, po > 0, but in all cases of the diverging 
system, eventually, p — > — oo, which is extreme non-SR behavior. 

In figure 6, another strongly divergent solution of the vector system is shown 
(see (69)). In spite of the divergence, the 3- velocity v x plotted seems to strongly 
converge towards a finite value. 

And finally, in figure 7, a case for an almost constant e = 0.5 case is shown. 
Eventually, the system decays due to the finite precision; it is therefore unstable. 

4.3 Analysis 

The converging cases depicted in figures (2,3,4) show, eventually, smooth 
uniform on-mass-shell motion. In particular, p, which is, essentially, the 3- 
acceleration squared af { as seen in a co-moving frame, is always positive, and 
eventually, falls to p — > 0, indicating uniform motion. 

The diverging case e — > oo requires a more detailed investigation, e 
indicates that \i\ 3> and therefore it follows that i 2 — x 2 = m 2 cg /M 2 cg >• 1. 

However 

e = i 2 -x 2 - l = i 2 (l-v 2 ) -1»0 (74) 

and therefore, v must be less than unity. 

At the present we do not have a way to deduce the actual limiting value 
of v in a diverging scenario, as the numerical simulation cannot continue the 
integration of a system where e — > oo indefinitely, though it does indicate a 
decreasing bound on it. 

It seems then, that even in the diverging e — > oo case, the 3-velocity con- 
verges towards a finite value, which would manifest itself as uniform motion in 
any frame, i.e., similar to the e — > case. 

Thus, the theory is well defined even for extreme unstable off-shell deviation. 

4.4 Origin of divergence 

The set of scalar equations (71) seems as a linear combination of the K t func- 
tions, coupled linearly to the dynamic variables. Therefore, it seems useful to 
study the behavior of the Ki functions, as possible origin of instability. Log- 
arithmically colored plots of the positive parts of Ki functions are shown in 
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Figure 5: A diverging solution example e — > oo, solved using the scalar system 
of equations (71). Initial conditions are £o = 0.5, po = —0.1 and all the rest are 
0. It strongly diverges at a finite r w 0.787. Only e, ln(e) and In \p\ are shown, 
though the behavior of the other terms are very similar, e, e and e are positive, 
whereas p, p, r\ are negative. 
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Figure 6: A diverging solution example e — > oo, formulated under the vector 
equation system given in (69). Only the plots for e, i, x and v x = x/i are shown. 
Interestingly, even though x and t diverge, the 3-velocity v x converges to a lower 
and finite value (less than unity). 
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Figure 7: A case of initially constant e in an accelerated system p = 0.1, 770 ~ 
—0.685, where, eventually, e — > possibly due to the finite, though large, nu- 
merical precision. The constant non-zero value of e appears to be unstable. 



25 



figures 8, 9 and 10, where we denned the notation: 

Kt=0{K i )K i = { Ki (75) 
V ' \0 Ki<0 

Of these functions, only Ki has substantial positive value even near the 
(e, e) — > (0, 0) origin' , as both K\ and K 3 have a factor of e. This can be seen 
as K 2 value depends linearly on the value of 4e + 3p, as well, and thus, even 
when i = 0, it may possess a positive value. 

The positive values of Ki normally push the system towards e > 0, i.e., 
towards the diverging e — > oo. 

To see this, we may begin with the equation for e in (71). The diverging case 
is characterized by e — > oo, and therefore, e, e, e -> oo as well. In the equation 
for e, we find that this divergence is indeed maintained, as long as p — > — oo, 
as well as iQ — > oo, and finally, the coefficient of K% has to be positive, leading 
to \p\ < e/2. I.e., though in this case p —> — oo, its divergence is controlled by 
p > -e/2. 

This characterization of divergence is also consistent with the equations of 
motion for p and 77, as long as r\ — > —00 for the diverging case. 
Summarizing, the conditions of instability are: 

Ki -> +00 

£,£,£•' ^+00 (76) 
p,p.T] -OO 

In figure 11, the value of the potential scalars Ki is shown for a converging 
case e — > 0. As r increases, eventually, i^i — >• 0,^2,^3 — > —00, or, K, t < 0. 
On the other hand, figure 12, a diverging system is shown where all the scalars 
Ki — > 00, and the entire system is attracted towards instability. 

4.5 Stability analysis 

Another way to see the divergence/convergence, is to use linearized stability 
analysis. The scalar equation set (71) can be cast into a formal first order 
system 

x = f(x) (77) 

where x € K 6 : 

x T = (e, i, e, p, p, 77) (78) 

7 The exact point e = is outside the domain of the above-shell analysis presented in this 
paper. 
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Figure 8: Logarithmically colored elevated plot of K*(e, e) (labeled as e and e' 
in this Visit [(>) plot), where Kf refers to the positive part of K^, the negative 
part being zero. 
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Figure 9: Logarithmically colored plot of K^{£, i, £, •) on a rectangular mesh 
(e, e and e are labelled as e, e' and e" in this Visit [()] plot). The valued nodes 
of the mesh were excluded to reveal the positive ones. Note that though the 
plot is 3D and shows only the dependence on e, e and e, the dependence on the 
extra p is co-linear with e, and therefore, there's no loss of information in this 
reduced function plot. The origin (e, e,e) = (0,0,0) is marked in the middle of 
the plot. 




Figure 10: Logarithmically colored elevated plot of K^(e, i) (labeled as e and e' 
in this Visit [6] plot). Note that the region e < is entirely zero, as K 3 (e, e) < 
there. 
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Figure 11: A case of a converging scalar system (71), where the value of the 
scalar potentials Ki are also shown. Clearly, all potential scalars Ki < as 
r — > oo, which leads to system stability. 
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Figure 12: A case of a diverging scalar system (71), where the value of the 
scalar potentials K t are also shown. Clearly, all potential scalars —¥ oo as 
t — > 1.6, in which the system becomes unstable. 
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and f(x) is a 6- vector defined as: 



-3/o + 2(e + l)Ki + sK 2 + (e + 2p)K 3 



P 



\ 



2n - iKi + 
-{e + 2p)K 1 



2pK 2 + pK 3 
f pK 2 + 2nK 3 



(79) 



The Jacobian J(x) is defined as: 

ja 
J b 



df a (x) 
dx b 



(80) 



On each point x along the orbit, J a b determines the local behavior, in the sense 
that for each eigenvalue with its associated eigenvector Vi, determine the local 
behavior of the motion. 

In figure 13, the real part of the eigenvalues of (80) along the orbit of a 
converging system are shown. In figure 14, on the other hand, the same plots 
for the case of a diverging system are shown. 

This suggests that the origin of the 6D phase space e, i, . . . , n — > is a stable 
stationary point, as eigenvalues near it are almost all negative. On the other 
hand, stationary points not near the 6D origin, are unstable. 

To see this, plots of the eigenvalues for an initially constant e are shown in 
figures 15, 16 and 17. 

In figures 15 and 16, the motion is of an accelerated particle such that its e 
is constant (see (73)). Clearly, in both plots, one of the eigenvalues is positive. 
Limited precision of the numerical simulation eventually causes the system to 
shift away from the constant e, leading either to a diverging e — > oo as in 15 
or converging e — > as in 16. Interestingly, even for uniform motion where 
e > 17, one of the eigenvalues is positive, which suggests that even this case 
is unstable under small perturbations. 



4.6 Dependence on M reg 

Both equation sets, i.e., equations (69) and (71), have a single parameter, D, 
defined in (62), which scales linearly with the renormalized mass D. In figure 
18, 4 plots of e(r) are shown, each having a different value of D. 

Thus, it appears as D 3> 0, the system becomes more stable, in the sense 
that conditions which cause diverging e — > oo cases, may be tamed to produce 
converging e — ¥ by increasing D. 

At first sight, it may seems contradictory that a large enough D would 
cause a given set of initial conditions to shift from divergence to convergence, as 
larger D increases K 2 , which would attract the system towards divergence even 
stronger. However, the key seems to be in the equation for e in (71). There, we 
find, for the D coefficient part 

7/2 

e = D 6 — [ - 16e(e + 1) + 16el + . . . = -16Dee 5/2 + ... (81) 
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Figure 13: Plot of the 6 eigenvalues of J a b of (80) along the orbit of the scalar 
system (71), in a case of convergence towards zero mass-shell deviation. Aside 
from one eigenvalue, all tend towards stable negative values as r —¥ oo. Only 
the real part is shown. 
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Figure 14: Plot of the 6 eigenvalues of J a b of (80) along the orbit of a diverging 
scalar system (71). 4 eigenvalues go towards +00, one goes towards —00, and 
one towards 0. 
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Figure 15: Plot of the 6 eigenvalues of J% of (80), for a scalar system (71), 
initially set for a constant e, where, eventually, due to limited precision, the 
system moves away from the stationary point towards divergence. This is one 
demonstration where constant e > cases are unstable. 
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Figure 16: Plot of the 6 eigenvalues of J% of (80), for a scalar system (71), 
initially set for a constant e, where, in this case, the limited precision leads to a 
convergence towards e — > 0. This is yet another demonstration of the instability 
of the constant positive mass-shell deviation. 
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Figure 17: Plot of the 6 eigenvalues of J a b of (80) for a scalar system (71) 
describing a uniformly moving particle with a positive constant mass-shell de- 
viation. 
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Figure 18: 4 plots of s(t), a solution of the scalar system (71), all starting with 
the same initial conditions. However, each system has a different parameter D 
(62), which scales linearly with the renormalized mass M reg . For the lower 3 
values of D, the system diverged, and for D = 1.3, the system converged. 
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I.e., in this case, larger D would cause a stronger attraction towards negative e. 
and therefore, towards lowering e. 

5 Conclusions and further study 

We have developed the renormalized Lorentz force (69) of a back-reacting par- 
ticle, using the r-retarded Green- Functions (38), first introduced in [2]. The 
original Lorentz force (14) was shown to depend on the sign of the mass-shell- 
deviation e (46) at and near the present point of the particle x^{t). In this 
paper, we studied the case of e > 0, in which the particle had a history, pos- 
sibly short one, in which it had a positive mass-shell deviation. In this case, 
even though the support of the Green-Function (38) is the entire past 5D light- 
cone, the equation was renormalized to avoid infinite force, after which a local 
equation remained. 

Note, however, that renormalization is not the only method. Other authors 
have chosen a different subtraction scheme than Gel'fand (see [9] and [17]), in 
which the integration on the past remains. The method chosen here, however, 
directly relies on Gal'fand [10] regularization method, and seems more natural, 
with the price of mass-renormalization (55). In this sense, this does not differ 
from other schemes, as even the well studied classical radiation-reaction in 4D 
requires mass renormalization if one includes only the retarded part of the 4D 
Maxwell green-Function (cf. [26, 23, 31]). The subsequent renormalized Lorentz 
force is 4 th -order vector equation (69), in which x v depends on lower deriva- 
tives, and scalars constructed out of them e, e etc. In the case of no external 
(renormalized) fields, as was studied in this paper, the vector equation has 12D 
phase space constructed of and x M . 

Contracting with (69) with x^,x M and i' M resulted in a reduced 6D system. 
This can be understood as follows: the original phase space essentially comprises 
three 4-vectors. We can construct six scalars out of these vectors: the size of 
each vector, e.g. x^x^, and the 3 hyperbolic angles between these vectors, e.g. 
x^x^. This is, essentially, the reduced scalar system (71). 

The first scalar, however, was e = — x^x 11 — 1, in which this particular form 
identifies it as a (renormalized) mass-shell deviation (to 2 — M 2 )/M 2 . 

Both the vector (69) and the scalar (71) equations where numerically inte- 
grated, showing that the system either converges towards the mass-shell e — > 0, 
or, with more slightly extreme initial conditions, e — > oo. In both cases, the 
3- velocity v = x/t seems to reach a finite value, which indicates that to an 
observer measuring x and t of such a particle would see, eventually, moving in 
uniform motion. The constant e case seems to be unstable, in the sense that 
small perturbations would cause e to change, either towards the converging or 
diverging cases. 

The converging e — > 0, is a classical description of a particle smoothly re- 
turning to its own mass-shell due to self-interaction. This remarkable result 
provides a mechanism for the mass stability of a charged particle in the Stueck- 
elberg framework. As shown in the previous section 4.6, if the parameter _D, 
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which scales linearly with the renormalized mass, increases sufficiently, even a 
diverging system can become convergent. 

We have seen that a divergence in e(t) is accompanied by the existence of 
a large positive eigenvalue (of the Jacobian matrix (80)) for local Lyapunov 
stability. A large positive eigenvalue can also have the effect of causing a rel- 
atively rapid transition from a seemingly diverging behavior to an evolution 
which eventually converges towards the mass shell e — > 0. 

As mentioned above, this paper focused on the above mass-shell e > case, 
which seems to be consistent, in the sense that the self-interaction system main- 
tains e > even in the converging cases. 

The limit point e = 0, however, is not handled by the above-mass-shell case. 
Therefore, in order to complete the study radiation-reaction in this framework, 
one would require a similar development of the Lorentz force (41) for the on- 
mass-shell e — and the below-mass-shell e < cases, and using external fields, 
describe how a particle could move in and out of these domains. This will be 
the subject of a future publication. 



A Canonical regularization of divergent integrals 

In this section we provide a short overview of the regularization method de- 
scribed in Gel'fand [ ]. 

The function x\ = 0(x)x x is non-zero for positive x, where 0{x) is the step- 
function. When acting on a smooth bounded function, 4>{x), 



(xX,4>{x)) = / x x <p{x) dx, 



which is well defined for 5R(A) > —1. On the other hand, the expression can be 
rewritten as 



{x\,4>{x)) = / x 



3=0 J ' 



dx 



(x) dx 



b x+ l +1 + 



(82) 



where the right-hand-side is well defined for {5R(A) > —to} H {A ^ —1, —2, }. 

This suggests that, as a generalized function, xh_ can be defined by its action 
on any smooth bounded function <p(x), as given by (82). The result is a function 
of A defined for all 5R(A) > —to except at A = —1, —2, ... — to + 1 where it has 

(pU) (o) 

simple poles with residues ; — . This suggests that xh_ itself is a generalized 

J- 

function with simple poles given by 



Res x. 



= (-1) 



IV. 
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Similarly, given 2 smooth functions, </>(x) and R(x), we are seeking a regu- 
larized solution for 

(i?; A (x),0(x)) = f*jj^te (83) 

where, a is defined by R(a) = 0, and R(x) > for x G (a, b) 8 . Furthermore, the 
expansion of R{x) around x — a begins at some order to > 

°° (t — n\ n+m 

= E t 1 \\ A n( a )> (84) 
(n + m) 

where A n are possible functions of a, the lower bound (much like r in (41)), and 
where we take Aq > 0, in order for R(x) > near x = a. 
Let us then then define the remainder function T{x): 

(x-a)™ ^ (" + m ) ! (85) 

1 ^ / S ( X ~ a ) A , N 

to! (to + 1)! 

where T(a) > 0, and the upper bound of (83) is taken such that T(x) ^ 
throughout a < x < b. Similarly, the expansion for 4>(x) around x = a is 
assumed to start at some order I > 

oo , \n+l 

W = E f " + ni g "( fl ) ( 86 ) 

n=0 ^ '' 

and similarly, the remainder function 

*>*>=<^ = l B ^ + 'i^ B ^ + - < 87 > 

The integral (83) then becomes 



Defining 



we then have 



(x - a) mA • T A (x) ' 7 Q T A (x) (x - a) mA -' 



'In the meantime, we assume R(b) > 0. 
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In case mX — I € {1/2, 3/2, . . .} is non-integer, then the regularization given 
in (82) is well defined, where we take x+ — > (x — a)™ A ~' and (j>(x) — > q{x). 

When mX — I € N is an integer, then regularization no longer applies, and 
the result has a simple pole whose residue is given by 



Res (R+(x),(j)(x)) 



m\-l£N 



b\-iy 



V-I-l 

dx r ' 



(91) 



For the above- mass-shell deviation, e.g., as given in (54), we find A = 
5/2, m = 2 and I = 2, and therefore, mX — 1 = 3, which residue is the renormal- 
ized force given in (56). 
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